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Abstract 


25 At the land surface, higher soil moisture levels generally lead to both increased evaporation for a 

26 given amount of incoming radiation (increased “evaporation efficiency”) and increased runoff 

27 for a given amount of precipitation (increased “runoff efficiency”). Evaporation efficiency and 

28 runoff efficiency can thus be said to vary with each other, motivating the development of a 

29 unique hydroclimatic analysis framework. Using a simple water balance model fitted, in 

30 different experiments, with a wide variety of functional forms for evaporation and runoff 

31 efficiency, we transform net radiation and precipitation fields into fields of streamflow that can 

32 be directly evaluated against observations. The optimal combination of the functional forms - 

33 the combination that produces the most skillful streamflow simulations - provides an indication 

34 for how evaporation and runoff efficiencies vary with each other in nature, a relationship that can 

35 be said to define the overall character of land surface hydrological processes, at least to first 

36 order. The inferred optimal relationship is represented herein as a curve in “efficiency space” 

37 and should be valuable for the evaluation and development of GCM-based land surface models, 

38 which by this measure are often found to be suboptimal. 

39 

40 Capsule Summary: The relationship between evaporation and runoff at the land surface is 

4 1 captured in a simple “efficiency space” framework. 
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42 1. Introduction: Efficiency Space 

43 A land surface model, or LSM, is an essential component of a climate modeling system. 

44 The numerous studies that have quantified the impact of land processes on simulated climate 

45 (e.g., Shukla and Mintz 1982, Delworth and Manabe 1989, Dinneyer 2000, Koster et al. 2000, 

46 Seneviratne et al. 2006, to name a few; see Seneviratne et al. 2010 for a review) underscore the 

47 importance of realistic treatments of land processes in climate models. Recognition of this 

48 importance has spawned multi-national and multi-institutional LSM evaluation projects such as 

49 the Project for the Intercomparison of Land-surface Parameterization Schemes, or PILPS 

so (Henderson-Sellers et al. 1993, Chen et al. 1997, Wood et al 1998, Bowling et al. 2003), the 

51 Global Soil Wetness Project (Dinneyer et al. 1999, Boone et al. 2004, Dinneyer et al. 2006), and 

52 the more recent benchmarking project PLUMBER (PALS Land Surface Model Benchmarking 

53 Evaluation Project, Best et al. 2014). 

54 One of the findings from PILPS is that in terms of LSM behavior and performance, the 

55 evaporation and runoff formulations of an LSM are inextricably linked (Koster and Milly 1997). 

56 (Note that throughout this paper, the tenn ‘evaporation’ is used in place of the more cumbersome 

57 ‘evapotranspiration’ to encompass all evaporation processes.) A land surface scheme’s treatment 

58 of the evaporation process has a profound impact on the runoff it generates; similarly, the 

59 scheme’s treatment of the runoff process has a first-order impact on evaporation. Perhaps 

60 counter-intuitively, an LSM with an excellent evaporation fonnulation will, when forced with 

61 realistic meteorology, produce poor evaporation rates if the model’s runoff formulation is poor, 

62 and it will accordingly lead to biases in any atmospheric model that is connected to it. To 

63 produce a reliable treatment of land surface processes in an Earth system model, it is arguably 

64 critical to understand the evaporation-runoff connection. 
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65 The framework introduced herein can serve as a tool for investigating this connection. 

66 Building on the work of Koster and Milly (1997) and Koster and Mahanama (2012), the 

67 framework captures the essence of the evaporation-runoff connection in a very simple and 

68 concise way: as a curve lying in “efficiency space”. 

69 In essence, the efficiency space framework is built on the assumption that evaporation 

70 efficiency (the ratio of latent heat flux, AE, to net radiation, R ne t, where A. is the latent heat of 

71 vaporization) is a function of soil moisture, W: 

72 AE/R net =p(W), ( 1 ) 

73 and on the further assumption that runoff efficiency (the ratio of runoff production, Q, over some 

74 time period to precipitation, P) is also a function of soil moisture: 

75 Q / P = F (W). ( 2 ) 

76 Under these two assumptions, the two efficiencies (i.e., AE/R ne t and Q/P) can be related to each 

77 other, as illustrated in Figure la. The red curve on the right shows an arbitrary functional 

78 relationship between AE/R ne t and W, and the blue curve directly below it shows a similarly 

79 arbitrary relationship between Q/P and W. Note that for mathematical tractability, both 

so evaporation efficiency and runoff efficiency are assumed throughout this paper to increase 
si monotonically with soil moisture, reflecting the reasonable idea that wetter soils allow an easier 

82 production of both evaporation and runoff. To capture the well-known plateau for wetter soils 

83 associated with energy-limited evaporation (e.g., Manabe, 1969), the rate of increase at high soil 

84 moisture can be set very low, as in the figure. 

85 For a soil moisture of 0.3 (degree of saturation), these functions produce an evaporation 

86 efficiency of about 0.3 and a runoff efficiency of about 0.05. For a soil moisture of 0.5, the 

87 evaporation and runoff efficiencies are about 0.5 and 0.1, respectively. These pairings of 
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88 evaporation and runoff efficiencies, along with every other possible pairing determined from the 

89 two functions, are plotted in “efficiency space” on the left side of Figure la. The efficiency 

90 space plot essentially shows how evaporation and runoff efficiencies vary with each other - how 

91 efficient, for example, the land surface is at converting precipitation into runoff when it is able to 

92 convert, say, 40 % of net incoming radiative energy into evaporation, and how the runoff 

93 efficiency changes when the evaporation efficiency is reduced to 30 % through the drying of the 

94 soil. 

95 Again, the evaporation and runoff efficiency functions on the right in Figure la are 

96 arbitrary; a different set of functions would lead to a different curve in efficiency space. Figure 

97 lb shows another possible example. The number of possible efficiency function combinations, 

98 of course, is infinite. The corresponding Q/P-vs-kE/R net curves could lie anywhere in efficiency 

99 space, subject to the condition that Q/P and AE/R ne t increase monotonically with each other. 

100 This paper has two overarching goals: (i) demonstrating that the identification of a curve 

101 in efficiency space is tantamount to characterizing overall hydrological behavior, and (ii) finding 

102 the curve that best represents nature and thus best serves as a target for LSM development. It is 

103 important to emphasize here that simple functional forms such as those on the right in Figure la 

104 are, at best, crude first order approximations to the much more complex behavior seen in land 

105 surface models and (presumably) in nature itself. Nevertheless, the functional forms capture 

106 enough of the underlying controls of soil moisture on evaporation and runoff to make the curve 

107 on the left meaningful and to allow nature’s efficiency space curve, however it looks, to embody 

108 in a valuable and concise way an essential aspect of surface hydrology. 

109 Given the difficulty of measuring evaporation and soil moisture at the large scale, the 
no efficiency functions underlying nature’s efficiency space curve, and thus the curve itself, cannot 
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in be derived directly from observations. This paper thus employs an indirect approach to search 

112 for nature’s curve: the use of a simple water balance model in conjunction with streamflow, 

113 precipitation, and radiation observations. The approach and the main results are presented in 

114 Section 2. The paper continues in Section 3 with a discussion of the efficiency space framework 
ns in the context of Budyko’s (1974) landmark hydroclimatic analysis, and Section 4 discusses the 

116 use of the framework for LSM development. The summary section (Section 5) includes a 

117 discussion of the framework’s use in the context of the new soil moisture data expected from 
ns recent and upcoming satellite missions. 

119 Underlying all of the discussions in the paper is the fundamental idea that it is the joint 

120 response of evaporation and runoff production to variations in soil moisture that underlies 

121 hydrological behavior. Focusing inordinately on one or the other in LSM development is likely 

122 to lead to inadequate LSM performance. 

123 

124 2 . Analysis of Efficiency Space with a Simple Water Balance Model 

125 

126 a. Model Structure 

127 The tool used in this paper to evaluate different curves in efficiency space is the simple 

128 water balance model (WBM) of Koster and Mahanama (2012). The essence of the model is 

129 illustrated in Figure 2. In addition to a soil water holding capacity, the user chooses the model’s 

130 functional efficiency relationships [i.e., kE/R net = P(W); Q/P = F(W)]; shown in the figure are 

131 arbitrary sample functions. The model is then driven with observations-based daily precipitation 

132 and net radiation forcing over a long period of time, typically several decades. At a given time 

133 step, the WBM’s (3(W) function, in conjunction with the net radiation input and the current soil 
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134 moisture state, determines the evaporation for the time step. Similarly, the WBM’s F(W) 

135 function, in conjunction with the precipitation, determines the amount of runoff generated. The 

136 WBM’s soil moisture state is updated based on the prescribed precipitation rate and the 

137 computed evaporation and runoff values. Precipitation is added to a snowpack reservoir during 

138 subfreezing periods; snowpack water melts and is added to the soil (essentially as precipitation) 

139 as the prescribed temperature warms. 

140 This WBM is undeniably much simpler than the land surface models (LSMs) typically 

141 used in climate modeling systems. The WBM lacks temperature prognostic variables and 

142 associated energy balance calculations, and it includes no specific treatment of (for example) 

143 baseflow and interception loss. Also, unless otherwise stated, most applications of the WBM 

144 here utilize the same efficiency functions in every season of the year and across large continental 

145 domains (i.e., without regard to spatial variations in topography, soil texture, vegetation type, 

146 and other features that may affect efficiency - shown in some studies to be important [e.g., 

147 Yildiz and Barros [2007]). This simplicity, while inappropriate for an LSM working within an 

148 Earth system model, is nevertheless justified for the WBM given that it is designed to contain 

149 only the most important, first-order controls on evaporation and runoff production contained 
iso within a more complex LSM (generally implicitly in the latter - the net result of complex 

151 interacting parameterizations), thereby allowing the isolated study of these controls. The 

152 numerous features not included in the WBM are thus assumed, when acting in a full LSM, to 

153 induce mainly second-order modifications to the fluxes. While the appropriateness of such an 

154 assumption can be argued, and while more complex versions of the WBM can easily be 

155 envisioned, it must be remembered that each addition of complexity brings the WBM closer to a 

156 full LSM and thereby makes the WBM results more difficult to interpret. The WBM, despite its 
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simplicity, has been found in tests to reproduce successfully the first-order behavior of complex, 
state-of-the-art LSMs. (See Koster and Mahanama (2012) for one example.) 


b. Simulation Procedure 

Following precisely the procedure of Koster and Mahanama (2012), the WBM is run on a 
2.5° x 2.5° grid across the conterminous United States (CONUS) over the period 1948-2000. The 
interannually-varying precipitation rates used to force the WBM are from the dataset of 
Andreidis et al. (2005), and the net radiation forcing used consists of climatological seasonal 
cycles derived from the Surface Radiation Budget (SRB) dataset (Gupta et al. 2006). The final 
52 years of simulated streamflow data are evaluated against observed streamflow data, which 
consist of naturalized streamgauge measurements in several large hydrological basins (see 
Mahanama et al. [2012] for details). 

Figure 3a illustrates a sample simulation. The red and blue curves in the left panel show, 
respectively, the evaporation efficiency and runoff efficiency functions prescribed in the WBM. 
The WBM is driven across CONUS for 53 years, and the runoff rates produced in the grid cells 
lying within the Upper Mississippi Basin (upstream of the Grafton stream gauge site) are 
combined into a basin-average annual rate for each simulation year. These simulated annual 
rates are compared to observed streamflows in the rightmost panel of Figure 3a. This particular 
combination of efficiency functions is seen to lead to significantly underestimated streamflow. 

Figure 3b shows a second sample simulation. The efficiency functions applied here are 
clearly different, particularly that for evaporation efficiency - in this second simulation, the 
ability of the WBM to convert net radiation energy to evaporation is, in general, reduced. As a 
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179 result, the simulated streamflows for the Upper Mississippi have increased and are now, in fact, 
iso too high. 

181 

182 c. Collapse of Functions onto Efficiency Space 

183 Figure 4 illustrates a finding that underlies much of the analysis in this paper. Three 

184 pairings of efficiency functions are provided in the top row of the figure. At first glance the 

185 pairings look fundamentally different, and yet a closer look shows that all three map onto the 

186 same single curve in efficiency space (lower left panel). Furthermore, all three pairings, when 

187 implemented into the WBM, produce essentially the same hydrological behavior - the lower 

188 right panel of Figure 4 shows that the time series of runoff for the three cases is essentially the 

189 same. This example is, in fact, representative. In general, pairings of efficiency functions that 

190 produce the same curve in efficiency space generate essentially the same hydrological fluxes in 

191 the WBM. 

192 This result can be considered in the context of the water balance equation: 

193 P = E + Q + C w AW/At, (3) 

194 where C w is the water holding capacity of the WBM and At is the time period considered. Using 

195 ( 1 ) and ( 2 ), this can be written as 

196 P = E + P G[kE/R net ] + C w AW/At, (4) 

197 where G is the functional relationship between P(W) and F(W) in efficiency space. Now 

198 consider the idealized limit of zero water holding capacity. In this case, given that P and R ne t are 

199 inputs, the equation collapses into one that defines a unique value of E for a time step and, 
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200 accordingly, also a unique value for Q. In other words, regardless of the evaporation function 

201 employed - regardless, for example, of which red curve in the top row of Figure 4 is used - a 

202 given curve in efficiency space (i.e., a given G[kE/R net ] function) implies a single set of values 

203 for evaporation and runoff and thus the same overall hydrological behavior. 

204 Naturally, the water holding capacity tenn is not necessarily negligible, and this explains 

205 the slight differences seen in the simulated streamflow curves in Figure 4 . Again, though, the 

206 overall strong agreement between the streamflow curves in this example is found to be standard 

207 through a comprehensive series of tests - for practical purposes, pairings of functions that map to 

208 the same efficiency space curve are found to be hydrologically very similar, supporting the use 

209 of efficiency space curves as an efficient means of characterizing the joint behavior of 

210 evaporation and runoff functions. 

211 

212 d. Efficiency Space Curves 

213 i. Overall ranking. The efficiency space curve in the lower left panel of Figure 4 , regardless of 

214 the pairing of efficiency functions that produced it, leads to streamflow values that are too low 

215 relative to observations (lower right panel). This efficiency space curve is reproduced in Figure 

216 5 a as the yellow curve, with the yellow color chosen to reflect this level of streamflow 

217 simulation skill. (Although no actual skill numbers are shown here, as this discussion focuses 

218 only on the relative skill levels associated with different curves, the skill is quantified with a 

219 root-mean-square error, or RMSE, calculation against observed ratios of annual streamflow to 

220 annual precipitation in the basin.) A pairing of efficiency functions that leads to the blue curve 

221 in Figure 5 a produces WBM streamflow values that are even more inaccurate (this is indeed the 
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222 pairing underlying Figure 3 b), whereas a pairing that leads to the red curve in Figure 5 a produces 

223 a relatively realistic simulation. 

224 Using this color-coding scheme for streamflow simulation accuracy, and plotting the 

225 higher scoring curves on top of the lower scoring ones, Figure 5 b shows the results of testing 

226 more than 23000 curves in efficiency space. Figure 5 b (and corresponding plots for other basins) 

227 holds a substantial amount of infonnation and is indeed a centerpiece result of this paper. It 

228 shows, for example, that any curve passing through the northwest quadrant of efficiency space 

229 will produce a poor WBM simulation of streamflow. Most importantly, given that observations 

230 are used to evaluate the skill levels, it provides an indication of the curve that effectively 

231 operates in nature across the Upper Mississippi Basin. The darkest red curves correspond to the 

232 highest skill values; any dark red curve will produce an accurate simulation and is thus a 

233 reasonable estimate of nature’s curve. 

234 The highest scoring curve (again, not by a large margin) is highlighted with a heavy black 

235 line. Assuming that this particular curve reflects nature, we can infer the following: during dry 

236 conditions, runoff production does not begin significantly until the evaporation efficiency is 

237 about 0 . 3 , and the runoff ratio is not especially sensitive to changes in soil moisture (i.e., relative 

238 to the sensitivity of evaporation efficiency to soil moisture) until the soil is wet enough to 

239 support an evaporation efficiency of about 0.6, at which point the runoff ratio becomes strongly 

240 sensitive to soil moisture. Evaporation efficiency, in contrast, is strongly sensitive to soil 

241 moisture variations during dry conditions and is much less sensitive to soil moisture variations 

242 during wet conditions. 

243 
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244 ii. Bearing of curve position on simulated streamflows. Figure 6 is provided to illustrate how 

245 different curves in efficiency space translate to different simulations of streamflow. When the 

246 WBM is fitted with evaporation and runoff efficiency functions corresponding to the black curve 

247 in Figure 5 b, the resulting simulation of streamflow is, as expected, quite accurate (Figure 6 a). 

248 Again, though, any curve lying in the dark red region of the plot will produce a reasonably 

249 accurate simulation. The right panel of Figure 6 b, for example, shows the WBM streamflows 

250 associated with the white curve in the left panel of Figure 6 b, a curve which lies within the dark 

251 red region. The simulated streamflows are only slightly less accurate than those in Figure 6 a, 

252 and given presumed uncertainties in the streamflow observations and in the meteorological 

253 forcing data used, one simulation cannot be said to be truly superior to the other. 

254 Worth mentioning here is that the skill levels shown in Figure 5 b are robust. If only the 

255 first half of the meteorological forcing and streamflow observations are used to establish the skill 

256 scores, the resulting optimal curves (not shown) are found to be very similar to those in the 

257 figure and result in similarly accurate streamflow simulations during the second half of the study 

258 period. 

259 Figure 6 c shows that when the efficiency space curve passes through the northwest 

260 quadrant, the resulting streamflow simulations are grossly overestimated. It seems safe to 

261 conclude that this curve does not come close to representing nature. On the other hand, when the 

262 efficiency space curve lies far below the black curve, as in Figure 6 d, the simulated streamflows 

263 are significantly underestimated. Figure 6 e shows the results obtained for an insensitivity of 

264 runoff efficiency to soil moisture. The simulated streamflows do vary slightly from year to year, 

265 but as expected, their interannual variability appears significantly underestimated relative to 

266 observations. 
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e. Water Holding Capacity 

The two main user-defined elements of the WBM are the evaporation efficiency and 
runoff efficiency functions, which, as shown above, can be described together as a single curve 
in efficiency space. The remaining user-prescribed WBM parameter is the water holding 
capacity, C w . In the simulations discussed above, C w is set to 169 mm, based on a chosen depth 
of 0.5 m, an assumed wilting point at a degree of saturation of 0.25, and a porosity of 0.45. For 
completeness, it is worth examining how the WBM results vary with C w . 

The efficiency space skill analysis in Figure 5b is repeated in Figure 7 for depths of 0.1 
m, 0.5 m, 1. m, and 1.667 m; again assuming a porosity of 0.45 and a wilting point of 0.25, these 
depths correspond to C w values of 34 mm, 169 mm, 338 mm, and 563 mm. Here the skill is 
evaluated over all of the shaded basins in Figure 7a collectively, with the skill contribution from 
each basin weighted by basin area. 

Two features from Figures 7b-7e stand out. First, while the optimal (black) curves differ 
somewhat between the four panels, the regions of high skill within efficiency space (the darker 
red areas) are basically the same - the highest scoring curves tend to hug the x-axis for lower 
values ofkE/R ne t and then climb sharply into the plot’s northeast quadrant for higher values of 
kE/R net . This basic curve shape is also seen for individual basins and other depths (not shown). 
The estimation of a valid relationship between evaporation efficiency and runoff efficiency 
functions is largely insensitive to the choice of C w . 

The second interesting feature of Figures 7b-7e is the fact that higher skill levels are 
possible with a depth of 0.5 m or 1 .0 m than with a depth of 0. 1 m or 1 .667 m. (The four panels 
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289 use the same color bar to illustrate skill, with levels that slightly differ from those used in Figure 

290 5.) This result could be interpreted in terms of the depth of soil “actively” participating in the 

291 land-surface hydrological cycle - perhaps, in a gross sense, changes in the soil water storage 

292 within the top 0.5 m or 1 m of soil are most relevant to the interannual variability of surface 

293 hydrological fluxes. Curiously, this depth range is roughly consistent with nominal estimates of 

294 vegetation rooting depth. Such agreement is perhaps not a coincidence. 

295 

296 f. Discussion 

297 The patterns shown in the efficiency space skill diagrams are entirely determined by the 

298 imposed observations-based precipitation and net radiation forcing and by the streamflow 

299 measurements used to evaluate the WBM output. Thus, the precipitation, net radiation, and 

300 streamflow observations, when examined jointly with the WBM, can be said to contain 

301 important and otherwise hidden information: (i) the fact that (based on the shape of the high 

302 scoring curves) evaporation’s sensitivity to soil moisture variation is greater than that of runoff in 

303 dry situations, whereas the reverse is true in wet situations; and (ii) the soil depth of greatest 

304 relevance to the generation of hydrological fluxes is perhaps of the order of 0.5 m to 1 m (Figure 

305 7). Most importantly, the high-scoring efficiency space curves suggest how evaporation and 

306 runoff efficiencies vary with each other in nature. This joint variation is critical to defining a 

307 basin’s hydrological behavior, at least in terms of how much runoff it produces each year; any 

308 other version of the joint variation (i.e., any curve in efficiency space lying outside the high- 

309 scoring region) produces inaccurate streamflows. 
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310 The high scoring curves for individual basins (not shown) are similar to first order, and 

311 yet they do nonetheless show some differences - there is, of course, some spatial heterogeneity 

312 in how nature behaves. Preliminary analyses of the WBM results indicate, for example, that the 

313 efficiency space curves for mountainous basins tend to lie a little above those for flatter basins in 

314 the dry regime, supporting the intuitive notion that mountainous areas produce runoff more 

315 easily than flat areas. Also, preliminary analyses suggest that flatter areas can produce higher 

316 skill scores with deeper soil depths, consistent with the fact that deeper depths-to-bedrock are 

317 found in flatter areas (e.g., as seen in STATSGO data [NRCS, 2012]). While a proper analysis 

3is of the spatial heterogeneity of basin behavior is beyond the scope of this paper, the potential for 

319 such heterogeneity must be kept in mind when interpreting the results presented herein. 

320 

321 3. The Budyko Perspective 

322 Budyko (1974), building on analyses by Schreiber (1904) and OTdekop (1911), 

323 pioneered the joint analysis of evaporation and runoff response to precipitation and net radiation 

324 forcing. Although Budyko’ s analyses are focused on climatic means - how, for example, annual 

325 mean evaporation varies with annual mean precipitation and net radiation - rather than on the 

326 short time scale (daily time steps or shorter) of relevance here, a comparison of Budyko’s results 

327 with those seen in the efficiency space diagrams is of interest. 

328 Budyko’s equation for evaporation is: 

329 E / P = [ (Rnet/Pk) tanh (Pk/R net ) ( 1 - cosh(R net /Pk) + sinh (R net /Pk) ] 1/2 . (5) 
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This function is plotted in Figure 8a. If we let the variable D represent the dryness index (i.e., D 
= R net / kP), then Budyko’s equation can be rewritten: 

kE/Rnet = f(D); (6) 

that is, the climatological evaporation efficiency is some complex function of dryness index. 
Now consider that for long-term averages, soil moisture variations in the water budget equation 
are irrelevant, so that the water balance can be written: 

P = Q + E. (7) 

Combining (6) and (7) then provides an equation for the climatological runoff efficiency: 

Q / P = 1 - D f(D). (8) 

Together, (6) and (8) allow the equivalent of Budyko’s climatological equation to be 
plotted in efficiency space. The idea is simple: a given value of D leads, through the two 
equations, to a unique pairing of kE/R net and Q/P and thus to a unique point in efficiency space. 
Examining a wide range of D values leads to a wide range of points, allowing the full curve to be 
plotted within the space. 

This curve is shown in white in Figure 8b. The color-coded skill levels appearing in the 
plot are taken from Figure 7c (i.e., they represent the WBM-based skill levels over the gray- 
colored basins in Figure 7a for a depth of 0.5 m). Two aspects of the Budyko curve in Figure 8b 
stand out. First, the curve has the same basic shape as the WBM-based high-scoring curves; 
indeed, Budyko’s curve falls within the high-scoring region of the space, though only just barely. 
Thus, even though Budyko’s equation is built to capture climatological efficiencies and 
implicitly (if empirically) accounts, for example, for the effects of seasonal phase differences in 
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351 the precipitation and net radiation forcing (Koster et al. 2006), the equation produces a curve in 

352 efficiency space that is not grossly inconsistent with the WBM-based curves, which focus on the 

353 short time-scale problem. 

354 Second, the Budyko-based curve is nevertheless seen to overestimate evaporation 

355 efficiency at the wet end, if we assume that either the optimal curve (in black) or the median 

356 curve within the high scoring region of the space is most correct. When the soil is fully wet (at 

357 the northeast comer of the plot), the Budyko curve implies that all of the net radiation will be 

358 converted into latent heat flux, i.e., that the land surface provides no resistance whatsoever to 

359 evaporation, even given the complexity of transpiration pathways through the vegetation. If such 

360 an extreme were true for climatological means, arguably it would also need to be true at any 

361 smaller time step, thus contradicting the median behavior of the WBM results. Some past 

362 studies (e.g., Koster et al. 2006) have indeed found observational evidence that Budyko’s 

363 estimates of climatological evaporation at the wet end are too high. 

364 

365 4. Relevance to Land Model Development 

366 

367 a. LSM Evaluation and the Model Dependence of Soil Moisture 

368 When evaluating the realism of an LSM, a key element to consider is the manner in 

369 which it allows variations of soil moisture to affect the magnitudes of evaporation and runoff 

370 fluxes. Figure 9 shows the diagnosed relationship between evaporation efficiency and soil 

371 moisture for three separate and representative LSMs (names withheld). The abscissa of each 
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point in each scatter plot corresponds to a monthly-averaged root zone soil moisture, and the 
ordinate of each point is the ratio of monthly-averaged evaporation to monthly-averaged net 
radiation; the data, taken from archived output of multidecadal offline simulations, are plotted 
here for points in and about the CONUS area during wann season months (May - September). 

A significant amount of scatter is seen in these plots, reflecting the fact (Section 1) that LSMs 
compute evaporation and runoff using a complex set of interacting parameterizations and thus do 
not generally use such efficiency functions directly. (These plots indeed demonstrate that the 
WBM is substantially simpler than a standard LSM.) Nevertheless, despite the scatter, a first 
order relationship between kE/R net and soil moisture can be seen for each model. 

Figure 9 shows why evaluating such relationships directly can be difficult. Within Model 
A, for example, the minimum soil moisture achieved varies spatially due to spatial variations in 
soil and vegetation properties, such as wilting point; for this model, there is no single relevant 
evaporation efficiency function. More importantly, evaluating an LSM’s efficiency function is 
difficult because the LSM’s soil moisture variable is model-specific, more of an “index of 
wetness” than a quantity that can be compared directly to observations (Koster et al. 2009). Such 
model dependency is reflected, for example, in the contrasting soil moisture ranges of Models B 
and C. Soil moisture and evaporation observations at the large scale are notoriously sparse, but 
even if they were plentiful, using them to evaluate a model’s evaporation efficiency relationship 
would be hampered substantially by this model dependency. 

The efficiency space framework, however, provides an elegant means for avoiding this 
problem. According to this framework, the model-dependent shapes of an LSM’s kE/R ne t-vs-W 
and Q/P-vs-W relationships are of secondary importance; the key relationship defining the 
LSM’s hydrological behavior is that between evaporation efficiency and runoff efficiency, a 
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395 relationship that can be examined without the direct consideration of soil moisture. The LSM 

396 can thus be evaluated by comparing its effective efficiency space curve against the 

397 corresponding observations-based curve (e.g., the black curve in Figure 7 c). 

398 

399 b. Using Efficiency Space to Evaluate State-of-the-Art Land Models 

400 Using archived model output to detennine directly the efficiency space curve that best 

401 characterizes a given LSM’s behavior comes with difficulties. Chief among these is the lack of 

402 an optimal averaging period over which to accumulate the evaporation, runoff, precipitation, and 

403 net radiation diagnostics. If the averaging period is too short (e.g., a week or so), then the runoff 

404 associated with baseflow, with its associated time delay, is no longer tied strongly to the 

405 concurrent precipitation. If, on the other hand, the averaging period is too long (e.g., a season), 

406 seasonal variability in soil moisture will confound the curve’s derivation - most of the runoff, for 

407 example, may be produced during a month of snowmelt, when soil moisture is high, whereas 

408 most of the evaporation may be produced a couple of months later, when the soil moisture is 

409 lower. Monthly diagnostics are a reasonable compromise, but even these, to some extent, are 

410 subject to such problems. 

4 11 To compute the LSM-based efficiency space curves shown in Figure 10 , a different 

412 approach is employed, one that avoids these problems and has the additional advantage of 

413 producing a simple, smooth curve in the space. In essence, the approach used to compute the 

4 1 4 highest scoring (black) curve in efficiency space in Figures 5 b and 7 is repeated, this time using 

415 the full LSM’s time series of streamflow values (rather than the observations) as the “truth” upon 

4 1 6 which to base the RMSE metric. The precipitation and net radiation forcing applied to the WBM 
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417 is derived from that used originally by the LSM in question, when it was generating its archived 

418 data. Thus, whereas the black curves in Figures 5b and 7c can be said to be in line with how 

419 nature works, the curve produced using the LSM-based “streamflow truth” can be said to be in 

420 line with how the LSM works - it can be said to characterize, in a gross sense, the LSM’s 

421 underlying relationship between evaporation and runoff efficiencies. 

422 The left panel of Figure 10a shows, in white, the optimal efficiency space curve 

423 established for a specific LSM (name withheld) based on its streamflow generation in the gray- 

424 colored basins of Figure 7a; plotted around it are the corresponding observations-based color- 

425 coded skill curves from Figure 7c. The curve for the LSM clearly lies above the high-scoring 

426 region of the space (where “nature’s curve” lies), especially toward the dry end. Based on the 

427 discussion of Figure 6c, this LSM can be expected to compute too much runoff. Sure enough, 

428 the right two panels in Figure 10a show that the annual runoff ratios computed by the full LSM, 

429 as detennined from the archived diagnostics, are too large compared to observations, with the 

430 largest errors seen toward the western half of CONUS. 

431 Note that a traditional interpretation of the comparison in the rightmost panels of Figure 

432 10a might be that the runoff formulation in this LSM is deficient. This may be true, but the 

433 efficiency space framework allows the deficiency to be stated in a potentially more helpful way: 

434 at the drier end, this LSM’s runoff efficiency is too large relative to its evaporation efficiency. 

435 Figure 10b shows the optimal curve obtained for a different LSM. Here the LSM-based 

436 curve falls below the high-scoring region, and, consistent with the discussion of Figure 6d, the 

437 resulting streamflows are underestimated relative to observations (rightmost panels). Again, the 

438 efficiency space framework allows us to describe this LSM’s deficiency in terms of its joint 
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439 evaporation and runoff behavior rather than solely in tenns of its runoff formulation. Correcting 

440 this LSM - bringing the white curve into the high-scoring region of the plot - could be a matter 

441 of improving its runoff fonnulation (increasing Q/P for a given soil moisture), but alternatively, 

442 and perhaps just as effectively, it could also result from improvements in the evaporation 

443 fonnulation (decreasing AE/R ne[ for a given soil moisture). 

444 Six other state-of-the-art LSMs were examined in the same way. The results (not shown) 

445 indicate that most of these LSMs are similar to Model Y in Figure 10b - most of the LSMs are 

446 characterized by efficiency space curves that fall below the high scoring region. 

447 

448 c. Potential for Simple Land Model Tuning 

449 Model Y in Figure 10b is in fact the Catchment-CN LSM (Koster et al. 2014), a merger 

450 of the energy and water balance framework of the Catchment LSM (Koster et al. 2000) and the 

451 dynamic phenology (prognostic biogeochemistry) components of the NCAR/DOE CLM4 

452 dynamic vegetation model (Oleson et al. 2010). As suggested by the figure, this relatively new 

453 LSM could benefit from further development or tuning of its evaporation and runoff 

454 fonnulations. 

455 When faced with the discrepancy in Figure 10b (i.e., the LSM’s optimal curve, in white, 

456 lying outside the plot’s high scoring region), one could envision two basic model development 

457 approaches. The first approach is simply to keep improving the evaporation and runoff 

458 fonnulations until the position of the white curve moves upward into the high scoring region. In 

459 principle this is the desirable approach, for if it is successful, the LSM’s curve would end up 

460 being in the right place for the right reasons. With this approach, knowledge of the relative 
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461 strengths of the evaporation and runoff fonnulations - knowing which formulation can be trusted 

462 more, based on the breadth of analysis underlying each - would point to the formulation that 

463 currently acts as the “weak link” in the simulation of hydrological fluxes and is thus the most 

464 appropriate target for focused development. The position of the white curve relative to the high 

465 scoring region would indicate a direction for the change in the model formulation. 

466 Improving a land model formulation through first principles, however, is far from trivial; 

467 depending on the accuracies desired, it may encompass an entire scientific career. Given time 

468 constraints and other considerations, the second model development approach - tuning an LSM 

469 fonnulation using the observations-based efficiency space curves - may, at least for some 

470 applications, be advantageous. 

471 The tuning approach in its simplest form is illustrated in Figure 11. The first step 

472 involves deciding whether the LSM’s evaporation or runoff fonnulation can be trusted the most. 

473 Compared to some LSMs, the Catchment-CN LSM has a relatively complex runoff calculation; 

474 still, the complexity built into the runoff fonnulation pales in comparison to that built into the 

475 evaporation fonnulation, and thus the latter is considered, for purposes of this demonstration, to 

476 be more trustworthy. Plotted in Figure 1 la are monthly (May through September) values of 

477 kE/Rnet versus monthly values of root zone soil moisture. The scatter here is similar to that 

478 shown in Figure 9; again, though, a first-order relationship can be seen between the two 

479 quantities. The red curve in Figure 1 la, obtained by averaging the LE/R ne t values over soil 

480 moisture bins, is used in this exercise to characterize the relationship. 

481 Figure 1 lb shows the observations-based efficiency space curves from Figure 7c. The 

482 black dotted curve is taken to be the “target” efficiency space curve for the tuned land model; 
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483 this curve, which falls wholly within the high scoring region, is used instead of the solid black 

484 curve because it accounts for a larger range of kE/R net values and thus works better for the tuning 

485 procedure. With the curves in Figures 1 la and 1 lb in place - the first showing a relationship 

486 between kE/R ne t and W, and the second showing one between kE/R ne t and Q/P - a relationship 

487 between Q/P and W is easily derived (Figure 1 lc). 

488 For the next step in the demonstration, the somewhat complex runoff formulation in the 

489 Catchment-CN LSM is completely stripped out of the model and replaced by a much simpler 

490 fonnulation, one in which the total runoff for a given time step is computed by multiplying the 

491 time step’s incident precipitation by a Q/P ratio taken from Figure 1 lc, based on the current 

492 value of root zone soil moisture. The results of this tuning exercise are shown in Figure 12. The 

493 tuning, in general, leads to improved streamflow simulation. The tuned LSM (Figure 12b) - 

494 which still, by the way, maintains most of its complexity (its diurnally-varying energy balance 

495 calculations, its spatial variation in vegetation, and so on) - produces annual runoff ratios that 

496 are, for the most part, more in line with the observations (Figure 12c) than those obtained with 

497 the original LSM (Figure 12a). 

498 This is, of course, a simple and rather heavy-handed sample demonstration. A more 

499 refined tuning exercise, one that isolates, for example, the particular facets of the runoff 

500 fonnulation that are known to be especially weak or arbitrary and then modifies only those 

501 facets, would in fact be much more satisfying and justifiable. Also, tuning the model separately 

502 for each basin would presumably produce better basin-by-basin results than those shown in 

503 Figure 12b. Given the difficulty of improving LSMs through first principles, a comprehensive 

504 tuning study may lead to a straightforward and more reasonable way of inducing an LSM to 

505 produce a more realistic simulation of hydrology. The present example serves merely to 
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506 illustrate how the information content of the observations-based efficiency space curves (e.g., 

507 Figure 7c) can work its way into an LSM fonnulation. 

508 

509 5. Summary and Discussion 

510 The typical land surface model, or LSM, is constructed as follows: the modeler produces 

511 the best evaporation formulation and the best runoff formulation possible and then combines 

512 them (along with formulations for various other land processes) into a single system. The 

513 modeler then tests this LSM against observations of (for example) surface turbulent fluxes and 

514 local streamflow rates and uses the test results, in conjunction with his or her knowledge of 

515 current model weaknesses, to adjust the individual LSM components. While on the surface this 

516 seems reasonable, the evaporation and runoff processes within an LSM are often developed and 

517 examined independently, thereby ignoring the fundamental truth emphasized in this paper: a 
sis first order understanding of the hydrological behavior simulated by an LSM requires the joint 

519 analysis of its evaporation and runoff formulations, so that the variations of evaporation and 

520 runoff efficiencies with respect to each other are, to first order, understood. 

521 To explore this joint variation, and to provide a potentially useful tool for developing, 

522 evaluating, and improving LSM formulations, the present paper introduces the concept of 

523 efficiency space. A simple water balance model (WBM) used in conjunction with precipitation, 

524 net radiation, and streamflow observations produces curves in efficiency space that illustrate, in a 

525 concise and simple way, how evaporation and runoff efficiencies tend to vary with each other in 

526 nature. Considering together all of the gray basins in Figure 7a, nature’s curve in efficiency 

527 space lies somewhere within the dark red region of one of the lower panels of Figure 7. In its 
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528 own way, such a curve is conceptually powerful, as it captures, to first order, the net effect of 

529 land surface processes on annual hydrological fluxes. 

530 Of course, such an analysis strategy does not come without its caveats. As discussed in 

531 the text, for example, the WBM’s structure is necessarily very simple, being designed to capture 

532 only the most important controls on the surface moisture fluxes. Another caveat, not yet 

533 discussed, involves the “control volume” assumed for the analysis. One holistic view of 

534 evaporation from the land surface interprets it as the net result of a variety of interacting 

535 processes in the soil, in the vegetation canopy, and in the atmospheric boundary layer, all 

536 evolving together (Betts 2004). Betts (2004) argues that addressing the evaporation problem 

537 through only a subset of these processes can have significant limitations. While this is certainly 

538 a valid and useful view, the analysis of evaporation in the present paper effectively centers on an 

539 alternative control volume, one that represents the soil and vegetation from the top of the canopy 

540 to some depth below the surface (say, a meter or two). In a sense, this control volume represents 

541 the operating world of the land surface model component of a climate model - a system that is 

542 forced from without by precipitation and incoming radiation and that must partition this 

543 incoming water and energy into a variety of fluxes and storage changes. The choice here of this 

544 alternative control volume, with its neglect of a detailed consideration of boundary layer and 

545 cloud processes, is guided - and is in fact made necessary - by the particular focus of this paper, 

546 namely, the examination of the dual role of soil moisture in controlling both evaporation and 

547 runoff (streamflow) production. 

548 Keeping these caveats in mind, the curve representing nature’s behavior in Figure 7c can 

549 be used as a target for LSM development. By noting where a model’s efficiency curve lies 

550 relative to that of nature, the modeler can adjust either the evaporation or the runoff formulation 
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551 - presumably whichever one is known to be weaker - in the right direction to bring the curves in 

552 line. The result should be an overall improvement in the simulation of streamflow and 

553 evaporation. Using the efficiency curve as a target allows the model developer to sidestep a 

554 commonly encountered problem in the evaluation of evaporation and runoff fonnulations: the 

555 fact that soil moisture in an LSM is more of a model-dependent index of wetness than a quantity 

556 that can be compared directly to observations. 

557 While focusing on efficiency space takes soil moisture out of the LSM development 

558 problem in a clean way, this must not be construed to imply that measurements of soil moisture 

559 are not useful for hydrological analysis and LSM development. The advent of satellite-based 

560 soil moisture measurements, particularly L-band estimates from the SMOS (Kerr et al. 2010) and 

561 SMAP (Entekhabi et al. 2010) sensors, opens the door to new and potentially powerful research 

562 into hydrological processes, including that associated with joint evaporation and runoff 

563 efficiencies. Consider, for example, the evaporation efficiency and runoff efficiency curves in 

564 Figure 13. The shapes of the functions here are arbitrarily drawn, but together they capture the 

565 basic shape of nature’s curve in efficiency space. As suggested in the figure, evaporation- 

566 induced changes in soil moisture during non-rainy periods should be smaller and more variable 

567 in time when the soil is in the drier regime, whereas rainfall-induced changes should be smaller 

568 and more variable in time in the wetter regime. One can imagine that by sorting soil moisture 

569 changes in the new satellite data into subsets associated with rainy and non-rainy periods and 

570 then examining these changes in the context of efficiency function pairings that correspond to 

571 nature’s efficiency space curve, one could establish where on that curve nature is operating at a 

572 given time - a non-traditional and potentially very useful way to characterize the current “overall 

573 hydrological state” of the measured region. Perhaps the analysis of the remotely sensed data 
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574 could even constrain better the curve’s location. Such research, of course, would have to address 

575 a number of issues, including the limited vertical penetration (relative to the root zone) of the L- 

576 band signal; still, the possibilities are intriguing. 

577 In closing, it is worth noting again that the joint consideration here of evaporation and 

578 runoff processes is in some ways not at all new, as it echoes, for example, the landmark work of 

579 Budyko (1971) and his predecessors. These earlier works focused on climatological means; the 

580 present paper, in a sense, discusses the ideas in the context of the short time scale and thus in the 

581 context of the operation of LSMs. The hope here is that the efficiency space framework, which 

582 succinctly captures the evaporation-runoff connection, can be used to promote LSM 

583 development. In the long run, such development should lead to improvements in overall Earth 

584 system model perfonnance. 

585 
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1. a. Illustration of how a representative evaporation efficiency function (XE/R net vs. W, the red 
curve on the right) and runoff efficiency function (Q/P vs. W, the blue curve on the right) 
combine to produce a single curve in efficiency space (the green curve on the left), b. A 
second sample illustration. 

2. Schematic of the water balance model (WBM) used in this study. The forcing variables, 
imposed daily, are precipitation (P) and net radiation (R ne t); the WBM uses the imposed 
evaporation and runoff relationships to compute the daily evaporation (E) and runoff (Q) as a 
function of its prognostic water content (W). (Reproduced from Koster and Mahanama 
[ 2012 ].) 

3. a. Illustration of a typical WBM simulation. Evaporation and runoff efficiency curves (the red 
and blue curves on the left, respectively) are implemented into the WBM, and the WBM is 
driven with gridded observational data over a number of decades. In the course of the 
simulation, the WBM produces gridded runoff values that are spatially aggregated to 
hydrological basins for comparison with observed basin streamflows (right panel), b. A 
second representative example. 

4. Demonstration that different pairings of evaporation and runoff efficiency functions lead to 
essentially equivalent streamflow simulations if their corresponding efficiency space curves 
are the same. The three pairings of efficiency functions in the top row (evaporation efficiency 
in red, runoff efficiency in blue) each map to the single efficiency space curve in the lower 
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left panel. WBM simulations with each pairing of functions correspondingly produce very 
similar simulations of streamflow (lower right panel). 

5. a. Three efficiency space curves, color-coded according to the degree to which corresponding 
WBM simulations reproduce observed streamflows in the Upper Mississippi Basin. The red 
curve corresponds to the most skillful simulation, and the blue curve to the least skillful 
simulation. Skill is quantified with an RMSE metric, though numbers are not shown in order 
to emphasize relative skill, b. Same, but with over 23000 color-coded curves overlain on the 
plot. The solid black curve is the curve associated with the greatest simulation skill. 

6. a. WBM simulation of streamflow in the Upper Mississippi (red curve in right panel) 
associated with the highest-scoring curve in efficiency space (black curve in left panel; see 
Figure 5b). The black curve in the right panel represents the observed streamflows for the 
basin, b-e. Same, but for the WBM simulation associated with the white curve in the left 
panel. 

Figure 7. a. Basins (in gray) examined collectively in the analysis. See Mahanama et al. (2012) 
for details on the basins and on the streamflow observations taken therein, b. Color-coded 
efficiency space curves as in Figure 5b, but for an assumed 0.1 m soil depth and for the 
collected set of basins. The black curve represents the curve associated with the highest 
simulation skill, c. Same, but for an assumed 0.5 m soil depth, c. Same, but for an assumed 
lm soil depth, d. Same, but for an assumed 1.667 m soil depth. 

8. a. Budyko’s semi-empirical relationship between the evaporation-precipitation ratio and the 
dryness index, defined as R net /Pk. b. Efficiency space curves from Figure 7c. A transform of 
Budyko’s relationship is shown in white. 
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9. Scatter plots of evaporation efficiency (ratio of monthly-averaged E to monthly-averaged R ne t) 
versus monthly-averaged root zone soil moisture, as derived from the May-September output 
diagnostics produced over CONUS by three state-of-the-art LSMs. 

10. a. Efficiency space curve (in white) for a specific LSM (“Model X”), overlain on the color- 
coded curves from Figure 7c. To the right are plots of annual runoff ratio across several 
basins as computed from Model X diagnostics and from observations, b. Same, but for a 
different LSM (“Model Y”). 

1 1. a. Scatter plot of evaporation efficiency (ratio of monthly-averaged E to monthly-averaged 
Rnet) versus monthly-averaged root zone soil moisture, as derived from the May-September 
output diagnostics produced over CONUS by the Catchment-CN LSM. b. Color-coded 
efficiency space curves from Figure 7c, with one of the curves highlighted as a heavy dotted 
line. c. Resulting relationship between runoff efficiency and root zone soil moisture. 

12. a. Annual runoff efficiencies (ratio of annual basin-averaged runoff to annual basin-averaged 
precipitation) for several CONUS basins, as produced by the Catchment-CN LSM in a 
multidecadal offline simulation, b. Same, but for a tuned version of the Catchment-CN LSM 
(see text), c. Same, but for observations. 

13. Sample evaporation efficiency and runoff efficiency functions, with highlighted regimes of 
interest. 
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736 Figure 1. a. Illustration of how a representative evaporation efficiency function (XE/R net vs. W, 

737 the red curve on the right) and runoff efficiency function (Q/P vs. W, the blue curve on the right) 

738 combine to produce a single curve in efficiency space (the green curve on the left), b. A second 

739 sample illustration. 
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Figure 2. Schematic of the water balance model (WBM) used in this study. The forcing 
variables, imposed daily, are precipitation (P) and net radiation (R ne t); the WBM uses the 
imposed evaporation and runoff relationships to compute the daily evaporation (E) and runoff 
(Q) as a function of its prognostic water content (W). (Reproduced from Koster and Mahanama 
[ 2012 ].) 
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757 


758 Figure 3. a. Illustration of a typical WBM simulation. Evaporation and runoff efficiency curves 

759 (the red and blue curves on the left, respectively) are implemented into the WBM, and the WBM 

760 is driven with gridded observational data over a number of decades. In the course of the 

761 simulation, the WBM produces gridded runoff values that are spatially aggregated to 

762 hydrological basins for comparison with observed basin streamflows (right panel), b. A second 

763 representative example. 
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769 

770 Figure 4. Demonstration that different pairings of evaporation and runoff efficiency functions 

771 lead to essentially equivalent streamflow simulations if their corresponding efficiency space 

772 curves are the same. The three pairings of efficiency functions in the top row (evaporation 

773 efficiency in red, runoff efficiency in blue) each map to the single efficiency space curve in the 

774 lower left panel. WBM simulations with each pairing of functions correspondingly produce very 

775 similar simulations of streamflow (lower right panel). 
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778 



779 Figure 5. a. Three efficiency space curves, color-coded according to the degree to which 

780 corresponding WBM simulations reproduce observed streamflows in the Upper Mississippi 

781 Basin. The red curve corresponds to the most skillful simulation, and the blue curve to the least 

782 skillful simulation. Skill is quantified with an RMSE metric, though numbers are not shown in 

783 order to emphasize relative skill, b. Same, but with over 23000 color-coded curves overlain on 

784 the plot. The solid black curve is the curve associated with the greatest simulation skill. 
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786 Figure 6. a. WBM simulation of streamflow in the Upper Mississippi (red curve in right panel) 

787 associated with the highest-scoring curve in efficiency space (black curve in left panel; see 

788 Figure 5b). The black curve in the right panel represents the observed streamflows for the basin. 

789 b-e. Same, but for the WBM simulation associated with the white curve in the left panel. 
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791 Figure 7. a. Basins (in gray) examined collectively in the analysis. See Mahanama et al. (2012) 

792 for details on the basins and on the streamflow observations taken therein, b. Color-coded 

793 efficiency space curves as in Figure 5b, but for an assumed 0.1m soil depth and for the collected 

794 set of basins. The black curve represents the curve associated with the highest simulation skill. 

795 c. Same, but for an assumed 0.5 m soil depth, c. Same, but for an assumed lm soil depth, d. 

796 Same, but for an assumed 1.667 m soil depth. 

797 


41 




a. 12 

1.0 

0.8 

£ 06 

0.4 

0.2 
0.0 

0 1 2 3 4 5 

Rne/ P X 

b. 

Runoff 
efficiency: 

Q/P 


798 

799 

800 Figure 8. a. Budyko’s semi-empirical relationship between the evaporation-precipitation ratio 
sol and the dryness index, defined as R ncl /PA. b. Efficiency space curves from Figure 7c. A 

802 transform of Budyko’s relationship is shown in white. 
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804 


805 Figure 9. Scatter plots of evaporation efficiency (ratio of monthly-averaged E to monthly- 

806 averaged R ne t) versus monthly-averaged root zone soil moisture, as derived from the May- 

807 September output diagnostics produced over CONUS by three state-of-the-art LSMs. 
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812 Figure 10. a. Efficiency space curve (in white) for a specific LSM (“Model X”), overlain on the 

813 color-coded curves from Figure 7c. To the right are plots of annual runoff ratio across several 

814 basins as computed from Model X diagnostics and from observations, b. Same, but for a 
sis different LSM (“Model Y”). 
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sis Figure 11. a. Scatter plot of evaporation efficiency (ratio of monthly-averaged E to monthly- 

819 averaged R ne t) versus monthly-averaged root zone soil moisture, as derived from the May- 

820 September output diagnostics produced over CONUS by the Catchment-CN LSM. b. Color- 

821 coded efficiency space curves from Figure 7c, with one of the curves highlighted as a heavy 

822 dotted line. c. Resulting relationship between runoff efficiency and root zone soil moisture. 
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823 

824 Figure 12 . a. Annual runoff efficiencies (ratio of annual basin- averaged runoff to annual basin- 

825 averaged precipitation) for several CONUS basins, as produced by the Catchment-CN LSM in a 

826 multidecadal offline simulation, b. Same, but for a tuned version of the Catchment-CN LSM 

827 (see text), c. Same, but for observations. 
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830 Figure 13. Sample evaporation efficiency and runoff efficiency functions, with highlighted 

831 regimes of interest. 
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